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Abstract. The prompt emission of gamma-ray bursts (hereafter GRBs) 
probably comes from a highly relativistic wind which converts its kinetic 
energy into radiation via the formation of shocks within the wind itself. 
Such "internal shocks" can occur if the wind is generated with a highly 
C*~) ' non uniform distribution of the Lorentz factor T. Taking into account 

such a variable distribution of T, we estimate the expected thermal emis- 
sion of the relativistic wind when it becomes transparent. We compare 
this emission (temporal profile + spectrum) to the emission produced by 
00 ' the internal shocks. In most cases we predict a rather bright thermal emis- 

sion that could easily be detected. This favors acceleration mechanisms 
for the wind where the main energy reservoir is under magnetic rather 
than thermal form. Such scenarios can produce thermal X-ray precursors 
^ ' comparable to those observed by GINGA and WATCH/GRAN AT. 
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1. The photosphere of a relativistic wind 

> ' 

, We do not discuss here the nature of the source responsible for the initial energy 

release leading to the GRB. We assume that a relativistic wind has emerged at 
radius r acc with a variable distribution of the Lorentz factor T(t) (with T ;> 100). 

The energy and mass injection rates are E(t) and M(t) = E(t)/T(t)c 2 . The wind 
production lasts for a total duration t w (total ejected mass : M w ). At r = r acc , 
the wind is still optically thick, whereas at the radius ns where the internal shock 
phase starts, it has to be already transparent, as indicated by the dominant non- 
thermal nature of the observed GRB spectra. The radius at which each layer 
of the relativistic wind becomes transparent can be estimated (see Daigne & 
Mochkovitch 2002) : photons emitted by a layer at r will successively cross all the 
layers emitted earlier by the source before they can escape from the relativistic 
wind at radius r csc (r) : as F > 1, we usually have r esc (r) S> r. The optical 
depth of a layer at radius r is approximatively given by r(r) — kE /87rc 3 r 2 r, 
where the opacity k is mainly due to the ambient electrons (Thomson opacity). 
Therefore the photospheric radius (r ~ 1) of the layer is given by : r p h ~ 
kE /87rc 3 r 2 . Fig. 1 shows this radius for an example where E = 10 52 erg.s -1 , 
t w = 10 s and the Lorentz factor is given in the same figure. 
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Figure 1. Left : initial distribution of the Lorentz factor used in our 
example. Right : the photospheric radius of each layer (solid line) and 
the corresponding radius at which photons escape (dotted line) as a 
function of a mass coordinate in the wind (solid line). 



2. Time profile and spectrum of the photospheric emission 

The photospheric luminosity and temperature of each layer can be estimated in 
the context of the fireball model (see e.g. Piran 1999) : 

L ph c E ~ 2/3 and kT ph ~ kT° , (1) 

where the initial fireball temperature is given by kT° = 1.3 K^Vi 1/2 MeV (m 
is the mass of the central black hole in units of 10 Mq). The emitted photons 
will closely follow a Planck distribution, so that the phostopheric emission leads 
to a spectrum and a count rate in an energy band [E±, E2] : 

dnP h (E) _ (l + z)L ph x x 2 ^ d cph _ L ph 1 /■*» x 2 dx 
dEdt 2"pianck(^T ph ) 2 expx-1 12 kT ph 2pi anck J Xl expx-l' 

(2) 

where x = (1 + z)E j kT ph , x lj2 = (1 + z)E lj2 / kT ph and X P i anc k = tt 4 /15. 



3. Comparison with the emission from internal shocks 

The internal shock luminosity is assumed to be Lis = f-y E where / 7 is the 
efficiency for the conversion of kinetic energy into 7-rays by internal shocks. 
This leads to a spectrum and a count rate in energy band [Ei, E 2 ] given by 

(3) 

where E p is the observed peak energy of the internal shock spectrum (#1,2 = 
Ei^/Ep) and B is the usual Band spectral shape (with / +o ° xB(x)dx = 2B an d), 



dEdt (l + z)E 2 l Band \E P J 1Z (l + z)E p J Ba nd J X1 



The expected photospheric emission of GRBs in the internal shock model. 
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Figure 2. Spectrum (left) and time profile (right) of the example 
presented in fig. 1. The time profiles are plotted in the 3.5-8.5 keV 
(top) and 50-300 keV (bottom) bands. The three curves correspond 
to the thermal photospheric emission (dashed line), the non-thermal 
internal shock emission (thin solid line) and the total (thick solid line) . 

characterized by a low and high-energy slope a and (3 (Band et al. 1993). 
Fig. 2 shows the spectrum and the time profile obtained in the example of fig. 1 
including both the internal shock and the photospheric emission (with a = —1.0, 
(3 = —2.25 and z = 1). The internal shock emission is computed with the model 
developped by Daigne & Mochkovitch (1998). The global observed peak energy 
of the non-thermal photons is 200 keV. 



4. Discussion 

The ratio of the photospheric over internal shock count rates is given by : 

r»ph /i i \ rn t P 2 x2dx 

r> _ G 12 _ i f- f -1 f-V 4 ,. 1 / 2 ! 1 + Z > tj P ^Band Jxi cxpx-1 , 

R12 ~CE /7ai 52 ^ 200 keV J Planck JJ B{x)d X ■ (4) 

This expression indicates that the photospheric thermal emission is usually dom- 
inant in the X- and 7-ray bands (see fig. 2) ! This is clearly in disagreement 
with the observations showing that the GRB spectrum is non-thermal in this 
spectral range. From eq. 4, there are two possible solutions to recover GRBs 
dominated by the non-thermal internal shock component : (i) increase of / 7 , 
i.e. more efficient internal shocks. However, our calculations show that even for 
/ 7 close to unity, the thermal component would still be easily detectable ; (ii) 
decrease of L p h and kT p h , i.e. a less hot and bright photosphere, compared 
to the prediction of the standard fireball model. If we define A as the fraction 
of the energy released by the source which is initially under thermal form, the 
photospheric luminosity and temperature are given by : 



^ X E (^y 2/3 and kTph ~ A 1 /* kT° ^) _2/3 , (5) 

V ?"acc / V r acc / 
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Figure 3. Same example as in fig. 2 but with a less hot and luminous 
photosphere computed with A = 0.01 (see text). 

so that the expression of R\2 in eq. 4 has now to be multiplied by a factor A 3 / 4 . 
In the standard fireball model A = 1 but much smaller values can be expected if 
a large fraction of the energy realeased by the source is initially under magnetic 
form. Fig. 3 shows exactly the same example as fig. 1, the only difference being 
that the photospheric emission has been computed with A = 0.01. A non- 
thermal burst is recovered in agreement with the observations. A small thermal 
precursor is still present in X-rays. It can be more intense for intermediate values 
of A, so that the photospheric emission predicted in this case could explain the 
X-ray precursors observed in a few GRBs. 

5. Conclusion 

We have computed the photospheric thermal emission of a relativistic wind 
with a highly variable initial distribution of Lorentz factor. We find that : (1) 
The photosphere predicted in the standard fireball model for GRBs is usually 
too hot and luminous. The thermal photospheric emission dominates the non- 
thermal internal shock emission, even for large internal shock efficiencies. This 
is in disagreement with the observations. (2) A less luminous photosphere is 
expected if the energy released by the source is initially injected mainly under 
magnetic rather than thermal form. In this case, we recover GRB spectra and 
temporal profiles dominated by the non-thermal internal shock emission. (3) We 
then obtain X-ray precursors which may correspond to the precursor activity 
observed in a few GRBs by GINGA and WATCH/GRAN AT. These resultats 
are presented in more details in Daigne & Mochkovitch (2002). 
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